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ABSTRACT 

Wc have explored the possibility of driving pulsation modes in models of sdO stars 
in which the effects of element diffusion, gravitational settling and radiative levitation 
have been neglected so that the distribution of iron-peak elements remains uniform 
throughout the evolution. The stability of these models was determined using a non- 
adiabatic oscillations code. We analysed 27 sdO models from 16 different evolutionary 
sequences and discovered the first ever sdO models capable of driving high-radial 
order ^modes. In one model, the driving is by a classical K-mechanism due to the 
opacity bump from iron-peak elements at temperature ~ 200,000 K. In a second 
model, the driving result from the combined action of K-mechanisms operating in three 
distinct regions of the star: (i) a carbon-oxygen partial ionization zone at temperature 
~ 2 10^ K, (ii) a deeper region at temperature ^ 2 10^ K. which we attribute to 
ionization of argon, and (iii) at the transition from radiative to conductive opacity in 
the core of the star. 

Key words: stars: oscillations - stars: variables: other 



1 INTRODUCTION 

Hot subluminous O-stars (hereafter sdOs) are evolved ob- 
jects on the way to becoming low-mass white dwarfs. Most 
sdOs are thought to have an inert carbon-oxygen core, with 
the luminosity provided by helium and hydrogen burning 
shells. The wide ranges of surface parameters, from 40 to 
100 kK for effective temperature (Toff) and 4.0 to 6.5 for 
surface gravity (log g) , spreads sdOs across a broad region in 
the HR diagram, that is crossed by post Asymptotic Giant 
Branch (post-AGB), post Red Giant Branch (post-RGB) 
and post Extended Horizontal Branch (post EHB) single 
star evolutionary tracks. This degeneracy in evolutionary 
paths plus possible binary star formation scenarios hampers 
reliable determination of sdO origins. This situation is not 
helped by the wide variety in sdO spectra: While they share 
a common characteristic of a He II 4686A line, a diversity of 
lines with different strengths of He I, He II and sometimes H 
are possible, rende r ing spectral classifi cation an open issue 
IJefferv et al.lll997l : IStroer et al.ll2007l ). 

The spectral analysis and determination of physical pa- 
rameters has the drawback of requiring the use of non- 
LTE atmospheric codes due to high effective temperatures, 
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whi ch yielded disagreements in the l arges t existent stud- 
ies (|Dreizler et al.l Il990l : iTheiU et ahl Il994l) un t il rec ently. 
Recent analyses of 13 sdOs by Stroer et al.l |2003) and 
iHirsch fc Hebed (|2008l ) showed that they segregate into two 
main groups regarding helium content: (1) helium- deficient 
sdOs with subsolar He abundances, no traces of carbon or 
nitrogen and a large spread in Tcff and log and (2) helium- 
enriched sdOs with supersolar He abundances and C and/or 
N lines present, spreading between ~40-50 kK in T^s and 
5.5 to 6 in log p. Group (1) sdOs are classified as post-EHB 
descendants of hot subdwarf B-stars (hereafter sdBs). To ex- 
plain the origins of group (2) sdOs, the authors invoke differ- 
ent channels: (i) me rger of two helium white dwarfs formed 
by binary evolution (|Webbinklll98i : iHan et al.ll2002l ). which 
explains the He enrichment, but does not quite reproduce 
the observed Teff-logg distribution, (ii) non-canonical evo- 
lution of a single star that experiences a phase of mix- 
ing from the H -shell into the envelope while on the RGB 
(|Sweigartlll997l ) and (iii) evolution of a RGB star whose en- 
velope is almost completely stripped by Roche- lobe overflow 
in a binary, pre-emptying the He - core flash (jPriebe et al.l 
Il999h . This last scenario has problems in explain ing the ob- 
served metal enrichment. (See lStrder" et al.l (|2007l ) and refer- 
ences therein for a thorough review of sdO formation mech- 
anisms). 
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Another way to attack the problem of untangling the 
different evolutionary scenarios proposed for sdOs, and in 
which this paper focuses, is the use of asteroseismological 
methods (provided, of course, that sdOs are pulsating stars). 
Asteroseismological analysis can yield information such as 
the mass of the star, envelope, and inner structure through 
the Brunt- Vaisala frequency and/or sound speed, as long as 
we can identify the oscillation modes present in the star. 

The first theoretical explorati on of the potentia l of sdO s 
as pulsators was carried out by ICharpinet et al.l (|l997al l. 
The stability analysis of full stellar sdO models correspond- 
ing to seven evolutionary sequences starting on the EHB 
found all sdO models st able. 

The previous year, ICharpinet et al.l (|l996l ) had theo- 
retically predicted the pulsator behaviour of sdBs due to a 
K-mechanism caused by iron enhancement in the envelope 
of the star. Since the observational confirm ation of the first 
pulsating sdB star bv lKilkennv et al] (1 19971), w e now know of 
about forty of these stars ( Oreiro et al. 20091 ) termed V361 
Hya, or more commonly EC 14026 stars, which coexist in the 
same area of the HR diagram with non-pulsating sdBs. They 
cluster in Teff at ~33 000 K and logg ~5.8 dex and show os- 
cillations with typical periods in the range 80-200 s, which 
have been identified with low-degree and low-radial order 
p-modes. Today, the sdB var iable stars also inc lude the long 
period pulsators V1093 Her (ICreen et al.ll2003f) also excited 
by the same k. mechani sm ([Fontaine et al.ll2003l l , a handful 
of hybrid pulsators (e.e:.[Oreiro et al.ll2005l ) and a unique He- 
sdB star (| Ahmad fc JeffervlbOOsi r For more details on sdB 
stars w e refer the reader to the excellent review bv l0stensenl 
(|2009l ). 

The goal of the present work is to com- 
plete the observational search for sdO pul- 
sa tors (iRodn'guez-Lopez. Ulla &: Garridd l2007l: 



Rodn'guez-Lopez et al.l l2005l : Rodriguez- Lopez et al.l 
20041 ) that proceeded in the period 2003-2006, and which 



has yielded non-conclusive results. The thrilling discov- 
ery of the first sdO pulsator, SDSS J160043.6-f 074802.9, 
(abbreviated J1600+0748) by a team of South African 
astronomers (|Woudt et al.l |2006| ). gave real support to 
our theoretical exercise. The recent discovery of p-mode 
oscillatio ns in sdO models tha t included radiative levitation 
of iron l| Fontaine et al.l |2008| ) is able to account for the 
the very fast multiperiodic oscillations observed in the 
range ~60 to ~120 s, hence establishing the same iron 
opacity mechanism found in sdB stars as responsible 
for sdO oscillations, and reinforcing the conclusion that 
radiative levitation and gravitatio nal settling are two key 
ingredients in driving pulsations dCharpinet et al.l Il997bl : 
ICharpinet. Fontaine fc Brassardll2009l ). 

Our aim is to investigate the possibility of driving pulsa- 
tion modes in theoretical models of sdO stars - coming from 
fully evolutionary sequences- with the aid of a non-adiabatic 
oscillations code. Thus, sdO models were specifically built 
using the evolutionary code jmstar (|Lawlor fc MacDonaldl 
I2OO6I and references therein). The code was used to further 
evolve these sdO models to more advanced state of evolu- 
tion, providing the tracks taken in the HR diagram. The sdO 
models are used as input to the GraCo non - adiabatic oscilla- 
tions code (|Mova. Garrido fc Dupretll2004l : lMova fc Garridd 
|2008| ). which predicts which oscillation modes of the models 
are excited. 
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Figure 1. Diamonds mark the location of the sdO models in the 
log g-Tefi plane. Filled diamonds correspond to models 8, 10 and 
15 which will be thoroughly described in following sections. The 
err or bars give the T^ff a nd logg determi nation for J 1600+0748 
by iFontaine et al ] l^i) (bottom) and iRodn'gue z-Lopez et al.l 
ll2009al) (top). A rough sketch of the sdO domain and the sdBV 
and GW Vir instability strips is also given. 



The paper is structured as follows: In Section 2 we give 
some general remarks about the evolution code, the struc- 
tural models and the non-adiabatic oscillations code. Section 
3 briefly describes the non-adiabatic stability analysis car- 
ried on the sdO models. Sections 4, 5 and 6 describe the anal- 
ysis for the most interesting models found. We summarise 
all the analysed models in Section 7, and the discussion and 
conclusions in Section 8. 



2 THE TOOLS 

2.1 The evolution code and the sdO models 

To have accurate theoretical models of sdO stars is of great 
importance in the subsequent pulsation analysis. The distri- 
bution of the thermodynamical quantities, such as pressure, 
temperature, density, and the chemical composition in the 
interior of the star, will determine the theoretical spectrum 
of modes computed by the non-adiabatic oscillations code. 

The evolution code jmstar was used to calculate a total 
of 16 initial sdO models corresponding to different evolution- 
ary sequences (Table [T] shows some of their relevant proper- 
ties). The models were built in a first step to account for the 
physical parameters of the pulsator sdOs candidates found 
in our simultaneous observational survey, and in a second 
step, trying to map out the sdO domain in the HR diagram. 
All of the models were created by evolving a 1 M0 star, 
starting on the pre-Main Sequence. The initial metallicity 
is solar (Z = 0.02) except models 7-9 for which Z = 0.05. 
The code follows the complete evolution of the star contin- 
uously to the sdO stage, and beyond. Fig. [T] roughly depicts 
the sdO domain in the HR diagram and the location of the 
constructed models. Their domain is located between the 
instability strips of the sdB variables and the GW Vir stars 
(low-degree g-mode pulsators excited by the C/O opacity 
mechanism) 
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Table 1 . Main physical parameters and mass fractions of the sdO models (models named . 1 means they belong to the same evolutionary 
sequence). X(other) refers to the mass fractions of all the other elements. Z is the current metallicity of the model (in brackets the initial 
metallicity, when different, at the ZAMS). Models between horizontal lines have been calculated with subtle differences in the code of 
evolution. (Model names are chosen mainly according to their mass loss rate parameter and effective temperature). 
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X(H) 


X(He3) 
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X(C) 


X(N) 


X(0) 
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Z 


& name 


(K) 




(Mo) 




















1 p675.8057 


79 000 


5.70 


0.478 


0.675 


0.21 


6.1E-05 


0.72 


3.4E-02 


l.OE-02 


5.2E-03 


1.7E-02 


0.07 (0.02) 


1.1 p675.9065 


90 000 


6.50 


0.478 


0.675 


0.20 


5.4E-05 


0.73 


3.5E-02 


l.OE-02 


5.1E-03 


1.7E-02 


0.07 (0.02) 


2 p650_7960 


79 000 


5.95 


0.491 


0.650 


0.68 


2.4E-04 


0.30 


3.0E-03 


1.7E-03 


l.OE-02 


5.3E-03 


0.02 


3 etap 685 


55 000 


5.89 


0.471 


0.685 


0.43 


1.4E-04 


0.55 


9.6E-04 


8.1E-03 


5.6E-03 


5.3E-03 


0.02 


4 etap 690 


55 000 


5.95 


0.471 


0.690 


0.33 


4.5E-05 


0.65 


6.3E-04 


l.OE-02 


4.0E-03 


5.3E-03 


0.02 


5 etap 695 


55 200 


5.98 


0.471 


0.695 


0.27 


2.8E-05 


0.71 


9.0E-04 


l.lE-02 


3.2E-03 


5.4E-03 


0.02 


6 etap 700 


55 000 


6.02 


0.470 


0.700 


0.18 


5.6E-05 


0.78 


1.5E-02 


7.8E-03 


2.6E-03 


9.8E-03 


0.04 (0.02) 



7 eta 600t45 


45 


000 


4.95 


0.497 


0.600 


0.59 


2.5E 


04 


0.36 


7.3E- 


03 


4.2E 


03 


2.4E 


02 


1.5E 


02 





05 


7.1 eta 600170 


70 


000 


5.95 


0.497 


0.600 


0.59 


1.9E 


04 


0.35 


7.3E- 


03 


4.3E 


03 


2.4E 


02 


1.5E 


02 





05 


8 eta 650145 


45 


000 


5.26 


0.469 


0.650 


0.59 


1.2E 


03 


0.35 


7.4E- 


03 


4.2E 


03 


2.4E 


02 


1.5E 


02 





05 


8.1 eta 650170 


70 


000 


6.22 


0.468 


0.650 


0.59 


LIE 


03 


0.36 


7.4E- 


03 


4.2E 


03 


2.4E 


02 


1.5E 


02 





05 


9 eta 700t45 


45 


000 


6.13 


0.454 


0.700 


1.6E-03 


8.2E 


07 


0.93 


I.IE 


02 


2.7E- 


02 


3.7E 


03 


2.7E 


02 


0.07 


(0.05) 


9.1 eta 700170 


70 


000 


6.34 


0.453 


0.700 


1.63E-03 


3.6E 


07 


0.93 


I.IE 


02 


2.7E 


02 


3.7E 


03 


2.7E 


02 


0.07 


(0.05) 


10 eta 675t45mixi 


45 


000 


4.18 


0.479 


0.675 


0.26 


5.2E 


05 


0.60 


6.0E 


02 


9.5E 


03 


5.4E 


02 


1.7E 


02 


0.14 


(0.02) 


10.1 15773t54g45 


54 


500 


4.54 


0.479 


0.675 


0.24 


4.3E 


05 


0.61 


6.3E 


02 


9.7E 


03 


5.6E 


02 


1.7E 


02 


0.15 


(0.02) 


10.2 15873t58g47 


58 


300 


4.68 


0.479 


0.675 


0.24 


4.1E 


05 


0.62 


6.4E 


02 


9.8E 


03 


5.7E 


02 


1.7E 


02 


0.14 


(0.02) 


10.3 16973t69g51 


69 


000 


5.14 


0.479 


0.675 


0.19 


2.4E 


05 


0.65 


7.1E 


02 


l.OE 


02 


6.2E 


02 


1.9E 


02 


0.16 


(0.02) 


10.4 17773t77g54 


77 


000 


5.39 


0.479 


0.675 


0.16 


1.4E 


05 


0.67 


7.7E 


02 


I.IE 


02 


6.7E 


02 


2.0E 


02 


0.17 


(0.02) 


10.5 18273t84g56 


84 


500 


5.63 


0.479 


0.675 


0.13 


9.0E 


06 


0.69 


8. IE 


02 


I.IE 


02 


7.0E 


02 


2. IE 


02 


0.18 


(0.02) 


10.6 18573t91g59 


90 


700 


5.93 


0.479 


0.675 


0.11 


5.9E 


06 


0.70 


8.4E 


02 


1.2E 


02 


7.3E 


02 


2. IE 


02 


0.19 


(0.02) 


11 eta 675t45mixnmi 


45 


000 


4.36 


0.480 





0.48 


3.0E 


04 


0.45 


2.8E 


02 


5.9E 


03 


2.9E 


02 


l.OE 


02 


0.07 


(0.02) 


12 eta 675t70mixli 


70 


000 


5.16 


0.479 


0.675 


0.19 


2.3E 


05 


0.65 


7.1E 


02 


I.IE 


02 


6.2E 


02 


1.9E 


02 


0.16 


(0.02) 


13 eta 675t70mix2i 


70 


000 


5.20 


0.479 


0.675 


0.16 


1.6E 


05 


0.66 


7.5E 


02 


I.IE 


02 


6.6E 


02 


2.0E 


02 


0.18 


(0.02) 


14 eta 675t70mixnmi 


70 


000 


5.43 


0.480 





0.48 


3.0E 


04 


0.45 


2.8E 


02 


5.9E 


03 


2.9E 


02 


l.OE 


02 


0.07 


(0.02) 


15 eta 675t90mixi 


90 


000 


5.86 


0.479 


0.675 


0.11 


6.2E 


06 


0.70 


8.4E 


02 


1.2E 


02 


7.3E 


02 


2. IE 


02 


0.19 


(0.02) 


15.1 18742t98g63 


97 


600 


6.28 


0.479 


0.675 


0.09 


3.4E 


06 


0.71 


8.8E 


02 


1.2E 


02 


7.6E 


02 


2.2E 


02 


0.20 


(0.02) 


16 eta 675t90mixnmi 


90 


000 


6.58 


0.480 





0.48 


3.0E 


04 


0.45 


2.8E 


02 


5.9E 


03 


2.9E 


02 


l.OE 


02 


0.07 


(0.02) 



To obtain sdO models, the mass loss rate on the RGB 
is enhanced by a factor ranging from 1.60 to 1.85 compared 
to canonical Reimers' mass loss rates. Most of the hydrogen- 
rich envelope is removed just after or even before the helium 
core flash. At the highest RGB rates, the model stars leave 
the RGB and evolve to higher temperatures at constant lu- 
minosity, avoiding the AGB. They then suffer a helium-core 
flash just before or slightly after they enter a white dwarf 
cooling phase. After the core flash, the star is an EHB/sdB 
star burning He quiescently in the core. After core He ex- 
haustion, the star enters the He shell burning sdO phase. 
The calculations are stopped when the model stars reach 
the (final) white dwarf phase. A born-again phase due to a 
very late thermal pulse occurs in some of the sequences. In 
this case, flash-driven convection carries nuclear processed 
material to the surface of the star, enriching the atmosphere 
with He, C, N and O. For particular models with sdO char- 
acteristics (see Table [!}, we have performed a stability anal- 
ysis. 

The evolutionary tracks on which models 1-2, 3-6, 7-9 
and 10-11 lie are depicted in Fig.[2l|3l|4]and[5] The complex 
evolution, near Toff = 30 000 K (log Toff ~ 4.5) in models 3 
to 6, and in a wider range of temperatures for models 7-9, 
occurs near the end of the helium-core flash when the star 
is settling to a quiescent helium burning phase in the core, 
that is, the sdB phase. At the high mass loss rates used to 



create models 3-6, the mass loss induced by the hot star 
wind is sufficient to expose internal He-rich layers. Thus, 
although the photospheres are not He-rich when the models 
reach the Horizontal Branch, they become He enriched due 
to the stellar wind as they evolve from sdBs to sdOs. Note 
that Fig. [2] shows only the part of the evolutionary tracks 
that occurs after the EHB phase. 

There is an important difference between the evolution- 
ary sequence of model 6 and the similar sequences of models 
3-5. For model 6, the He-core flash occurs after the star has 
entered the white dwarf cooling phase, and in a similar man- 
ner to the evolution of post- AGB 'Born- Again' stars (see e.g. 
iLawlor fc MacDonaldll2006l and references therein) , the con- 
vection zone that develops above the off-center He-burning 
region penetrates into the H-rich envelope. However, unlike 
the 'Born- Again' stars, this occurs when the He-burning re- 
gion has become inactive and the temperature is too low for 
proton capture reactions to occur. The result of the mixing 
is to introduce C^^ produced by the 3-a process into the 
envelope, together with O^* and Ne'^^ produced by N^* a- 
captures. Shortly after the end of the He-core flash, this ma- 
terial is convected throughout the H-rich envelope, increas- 
ing the heavy element abundance to Z = 0.04. For models 
3-5, the convection zone above the off-center He-burning 
region does not penetrate into H-rich envelope and there is 
no increase in the surface heavy element abundances. The 
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Figure 2. Evolutionary tracks from the horizontal branch de- 
scribed by models 1 and 2. A stability analysis was done also for 
a subsequent model in the evolutionary track of model 1, namely, 
model 1.1. Diamonds mark the location of the models. 
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Figure 4. Complete evolutionary tracks for models 7-9. A sta- 
bility analysis was done also for subsequent models in the evolu- 
tionary track, namely, models 7.1, 8.1 and 9.1. Diamonds mark 
the location of the six models. 
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Figure 3. Complete evolutionary tracks for models 3-6. Dia- 
monds mark the location of the models. 

temporal evolution of the photospheric abundances for the 
model 6 sequence is shown in Fig. |6] 

Models 10 to 16 require special comment. They are also 
late core helium flashers, but after they had arrived on the 
Zero Age Horizontal Branch (ZAHB) ad hoc extra-mixing 
was included for the purpose of bringing heavy elements to 
the surface to increase the envelope metallicity. Due to the 
extra-mixing, helium is also taken to higher temperatures, 
leading to helium shell flashes. Flash-driven convection then 
leads to further mixing of heavy elements to the surface. 
These helium shell flashes are responsible for the loops in 
Fig. El near T^s ~ 60 000 K and logg - 5. For models 11 
(see Fig. [SJ, 14 and 16 mass loss was switched off in the 
calculations once they arrived on the ZAHB. 

2.2 The oscillations code 

We used the GraCo fully non-adiabatic oscillations code 
(|Mova et al.ll2004l : lMova fc Garridoll2008r ). which is based on 




5.5 5.0 4.5 4.0 3.5 3.0 



log Teff (K) 

Figure 5. Evolutionary tracks from the horizontal branch for 
models 10 and 11. Dots mark the location of the models. 



the lUnno et~all (| 19891 ) formulation to study the sdO models 
described above. Because the thermal and dynamical time 
scales become comparable in the outer parts of the star, a 
non-adiabatic treatment of the oscillations is required in or- 
der to properly describe the dynamical characteristics of the 
oscillations, i.e. which modes will grow or decay with time. 
The non-adiabatic treatment does not signiflcantly modify 
the eigenfrequencies with respect to the adiabatic, at least 
for our purposes, although it will modify the amplitudes and 
phases of the eigenfunctions. This proves to be very impor- 
tant in using multicolor photometry for mode identiflcation 
of pulsating stars, although this is beyond the scope of this 
paper. 
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Figure 6. Temporal evolution of the photosphcric abundances 
of the model 6 sequence. Time is measured from the peak of the 
He-core flash. The vertical line indicates the location of the sdO 
model 6. 

3 NON ADIABATIC STABILITY ANALYSIS 

We used GraCo to calculate modes from I = to I = A 
with frequencies between ~0.3 to ~20 mHz (~50 to 3000 s), 
searching for excited short and/or long periods. 

The stability or instability of a model is given by the 
growth rate as a function of the frequency. We used the 
iDziembowski. Moskalik &i Pamvatnvkhl (|l993h definition of 
the growth rate: -q = W/ \dW/dr\dr with the convention 
that positive (negative) values indicate driving (damping) 
of a particular mode. 

In discussing which regions of the model contribute to 
driving or damping, we made use of the derivative of the 
work integral as a function of the fractional mass deptlQ 
(dW/dlogq), which gives the energy locally gained or lost by 
the displaced material during one pulsation cycle. A positive 
(negative) value of —dW/d\og g at a given location indicates 
that this region contributes locally to driving (damping) a 
mode. Correspondingly, the running work integral {—W) is 
also used to illustrate the relative importance of driving and 
damping zones. 

In the next sections we describe the more interesting 
models found in the analysis. 



4 MODEL 8: ETA 650T45 

This model, taken from an evolutionary sequence with ini- 
tial heavy element abundance Z = 0.05, has Toff =45 000 K 
and logg — 5.25. Although it was found stable in the 

^ The fractional mass depth parameter logg = log(l — Mt/Mt) 
gives a bettor resolution of the envelope structure, which is more 
relevant to pulsations than the structure of the deep interior. 



Model eta650t45 1 = 2 



. 'U ' ' 1 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ',. 

1 \ 1 

k^AAA A A. . I.\ ^s'. . ... ^ .... ^ A- A. A 


1 1 1 1 1 

.. . p30 - 
• • • - 




















I - \ \^-^ 

i 1 \ 




g88 - 


- 1 \ 

- 1 \ 

' > 







-2 -4 -5 -8 -10 -12 -14 

log q 

Figure 7. Propagation diagram: Brunt- Vaisala (solid line) and 
Lamb frequency (long-dashed line) for I = 2 modes of model 8. 
The opacity (dotted line) scale is given in the right axis. The 
vertical dashed lines mark the maximum of the C/O-He and He- 
ll chemical transitions. The frequencies of a pure g-modc (g88), 
a mixed mode (pi) and an almost pure p-mode (p30) are given, 
with the positions of their nodes marked by circles. See text for 
further details. 

whole frequency range, it has a narrow intermediate fre- 
quency region, corresponding to low-radial order p-modes, 
with a weak tendency to instability, which seemed interest- 
ing in terms of explaining a possible selection mechanism 
for observed modes in the only sdO pulsator found to date 
(|Rodrfguez- Lopez et all (|2009bl ). 

The model propagation diagram, representative of all 
the models used in the analyses, plus the opacity are shown 
in Fig. [7] There are three main bumps in the Rosseland 
opacity profile (scale given in the right axis): the shallowest 
at logq ~ — 12 corresponds to the helium ionization zone, 
where He II turns to He III at ~ 40 000 K; the maximum 
at logg ~ —9.5, which we will refer to as the Z-bump, is 
mainly due to iron-peak elements at ~ 200 000 K; and the 
one at log g ~ — 5 is the deep opacity bump at ~ 2 000 000 K, 
which results from a combination of K-shell photoionization 
of carbon and oxygen a nd spinorbit effects in L- shell bound- 
free transitions of iron (|Rogers fc Iglesiaslll99d ). It is clear 
from the opacity profile that the He II partial ionization zone 
is too high in the envelope to have any significant weight in 
the driving of the modes. 

The Brunt- Vaisala (logiV^) and the Lamb frequencies 
(logLf) are plotted for 1 = 2. The Brunt- Vaisala frequency 
has maxima where gradients in chemical composition are 
steepest. Hence, the peak at logg ~ —0.84 accounts for 
the transition between the C/O core and the He burning 
shell, while the secondary peak at log q ~ —2.23 is the tran- 
sition between the radiative He and H burning shell. All 
the energy transport in the model is radiative except at 
—9.4 ^ log q ^ —10.0 where A*'^ < 0, indicating that this is a 
convection zone associated with ionization of heavy elements 
producing an increase in the opacity rendering convective 
energy transport more efficient than radiation. We note here 
that models of sdB stars usually have a convective core and a 
radiative envelope with a thin conve ction zone produced b y 
the partial ionization of HcII/HcIII (jCharpinet et al.]|2000l ). 
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Figure 8. Left, Centre & Right: Differential work (solid line), running work integral (dotted line) and opacity (dashed line) for the glOO, 
gl5 and gl mode, respectively for a quadrupole mode of model 8. All plots have been scaled to arbitrary units. 




Figure 9. Left, Centre & Right: Differential work (solid line), running work integral (dotted line) and opacity (dashed line) for the p2, p9 
and pll mode, respectively for a quadrupole mode of model 8. All plots have been scaled to arbitrary units. The vertical dashcd-dotted 
lines in the last plot depict the convection zone. 
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Figure 10. Growth rate vs. frequency for 1 = 2 modes of model 8. 
Diamonds indicate the modes for which the differential work and 
running work integral were plotted. 

The g- and p-mode propagation zones are determined 
respectively by cr^ < N^,Lf and > N^,Lf, with a be- 
ing the real part of the eigenfrequency. Fig. [7] shows that 
g^modes cannot propagate in the outer envelope, and that 
p-modes cannot propagate in the innermost part of the star, 
except for modes with very high-radial ordeiQ. The diagram 
also gives broad constraints on the range of frequencies in 
which modes can propagate. It gives an upper limit for the 

^ The radial order, represented by k, or sometimes n, gives the 
number of nodes of each mode, in the radial direction. 



angular frequencies of modes, loga^ ^ ~1 ~ depending 
somewhat on the degree I - corresponding to cyclic frequen- 
cies < 50 mHz or periods P > 20 s. In the case of p-modes 
it gives a lower limit for eigenfrequencies, logo"^ > —4.3, 
which translates to cyclic frequencies i/ > 1.1 mHz or periods 
P < 890 s. Hence modes with oscillation eigenfrequencies in 
the range —4.3 < logcr^ < -1 are mixed modes, i.e. they be- 
have as p-modes in the outer part of the star and as g^modes 
in the inner zones. These limits correspond to a range in the 
frequency propagation of 1.1 < < 50 mHz or a range in 
periods of 20 < P < —890 s. These mixed modes are typical 
of evolved stars, like sdOs, and hamper the identification 
of the modes: P and G modes' propagation zones are not 
separated by an evanescent zone, but their frequency ranges 
overlap significantly, as derived above. We should then be 
cautious about the nomenclature and emphasise that when- 
ever we refer to a p- (g-) mode, we refer in fact to mixed 
modes with nodes mainly in the P (G) mode propagation 
zone (see position of nodes for modes pi, p30 in Fig. [7]); ex- 
cept for those with high-radial orders (see nodes for mode 
g88 in Fig. [7| i.e. g-modes with v < 1.1 niHz, and p-niodes 
with v ^ 50 mHz, far out of the frequency range considered 
here. In practice, for this model, modes with > 22 mHz 
(or k>33) do not have modes in the G mode region. 

We plot the growth rate rj vs. frequency for 1 = 2 modes 
(the behaviour being analogous for any degree I) in Fig. 1101 
up to 15 mHz, as higher frequencies were found highly sta- 
ble. The model is stable as the growth rate is negative in 
the whole frequency range. However there are 3 regions of 
interest in regard to stability: 

• between about 0.5 to 1 mHz, where modes are not com- 
pletely stable 
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• from about 2.5 to 5 mHz, a highly-stable region 

• from ~ 5 mHz to 15 mHz, where the growth rate shows 
an oscillatory behaviour. 

The differential work (— dVF/d log g), running work in- 
tegral { — W), and opacity, for six modes (marked with dia- 
monds in Fig. llUp chosen to be representative of the growth 
rate profile, give some insight into the local contributions 
and their relative importance to the work integral. 

Fig. |S] shows the glOO mode (on the left) with v = 
0.57 mHz and rj = —0.78, which is representative of the 
low frequency region where the growth rate is not so neg- 
ative. We identify a driving region at the location of the 
deep opacity bump as responsible for the not complete sta- 
bilization of the modes. However, due to the large number 
of nodes in the deepest region, radiative damping exceeds 
driving and the modes are stable. 

The gl5 mode (centre) with v = 2.91 mHz and 77 = —1 
represents the high-stability region found at frequencies be- 
tween 2.5-5 mHz. The region where maximum energy inter- 
change for the model occurs where the ratio of the dynamical 
to thermal timescales is ~ 1 (corresponding to logq ~ —10). 
For low-radial-order (?-modes, maximum energy interchange 
takes place at log q ~ —9, where the derivative of the opac- 
ity does not favour driving, yielding a complete stabilization 
of the modes. 

The onset of the oscillatory behaviour is represented by 
the gl mode (right) with v = 5.58 mHz and 77 = —0.81. The 
mode gains a bit in instability due to the development of a 
weak driving region at —9.5 < logg < —10, i.e. the Z-bump 
location. Its oscillating nature is due to the rise and fall of 
a driving zone at this Z-bump location, as we show next. 

We plot 3 modes belonging to the third region of sta- 
bility in Fig. [51 On the left panel, we plot p2 mode with 
V — 6.56 mHz and 77 = —0.98. We see that the reason for the 
low value of the growth rate is the vanishing driving region, 
whereas all the remaining energy contributes to damping. 

The mode with maximum growth rate is p9 (centre) 
with v = 9.30 mHz and rymax = —0.64. The reason is the 
maximum developing of the driving region at the location 
of the Z-bump. 

The decrease of the growth rate at high frequencies 
is due to the energy damping region becoming wider and 
wider, as is the case for the pll mode (right) with v = 
9.98 mHz and rj = -0.79. 

Thus, the interplay between the development and ex- 
tinction of a driving region at the Z-bump location is re- 
sponsible for the oscillatory behaviour of the growth rate. 

In an attempt to thoroughly investigate the properties 
of this model (and in particular the peculiar behaviour of 
the growth rate), we explored the behaviour of the kinetic 
energy of the modes. This lead to the discovery of a mode 
trapping structure in the g- and jp-m ode spectra, the details 
of whi ch will be presented elsewhere (jRodrfguez-Lopez et al.l 
l2009bl ). 



5 MODEL 10: ETA 675T45MIXI 

This is the first of seven models - models 10 to 16, see Ta- 
ble [1]- created by introducing an artificial extra-mixing in 
the post-ZAHB evolution with the purpose of raising the 



metallicity in the envelope in an attempt to drive pulsations, 
as explained in Section 2.1. 

Model 10 was our first model discovered to be unsta- 
ble. We found that high-radial order g^modes are driven 
through the action of a K-mechanism due to the Z-bump 
at ^ 200, 000 K, the same mechanism that drives oscilla- 
tions in sdB stars. Although the extra-mixing has signifi- 
cantly increased the abundances of C and O in the driving 
region, this has very little effect on the opacity there. Hence 
the presence of instability in this model shows that radia- 
tive levitation of iron is not a necessary condition for pulsa- 
tions of sdO stars, as long as the metallicity is sufficiently 
increased. More evolved models from this evolutionary se- 
quence (models 10.1 - 10.6 in Table [1} were also tested for 
pulsational instability. We found pulsations were driven in 
models with Tes ^ 54 000 K and logg 4.5. So, the insta- 
bility strip for this evolutionary sequence was found to be 
45 000 ^ T^ff < 54 000 K and 4.2 < logfl ^ 4.5, and may 
extend to lower temperature and surface gravity. 

The evolutionary track on the \ogg~TcS diagram de- 
scribed by the time evolution of the model is shown in 
Fig. 111! The positions of the initial model 10 and the sub- 
sequent model found unstable are marked with diamonds. 
Crosses indicate intermediate structural models found sta- 
ble in the stability analysis. The whole evolutionary time for 
this model is about 300 000 yr, while e-folding times of the 
excited modes are of the order of 10~^ to 10 yr, so modes 
have enough time to grow in amplitude. 

The growth rate profile for I = 2 modes is shown in 
Fig. 1121 We found two regions with excited modes (similar 
growth rate profiles and instability regions where found for 
the other degrees): 

• within the frequency range 0.29 ^ ^ 0.32 mHz (pe- 
riods 52-58 min) , which correspond to modes from g203 to 
gl82 and 

• within the frequency range 0.38 ^ ^ 0.42 mHz (pe- 
riods 40-44 min), which correspond to modes from gl52 to 
gl38 

These modes fall within a wider range favoured for in- 
stability between 0.2 < < 1 mHz. Meanwhile, modes 
with frequencies v > 2 mHz were found highly stable with 
ri = —1.0. We have selected a few representative modes, 
marked with diamonds in Fig. 1121 for further discussion. The 
profiles of the differential work, running work integral, and 
opacity for these representative modes are shown in Fig. 1131 

We chose the gl90 mode with 1/ — 0.31 mHz and 
ri = 0.17 (Fig. 1131 left) as representative of the two 'is- 
lands' of instability. There is a wide driving zone centered at 
logg — —8, the location of the Z-bump responsible for the 
instability. 

To represent the stability zone between the two 'islands' 
of instability and, moreover, the stability zone up to 2 mHz, 
we chose the gl77 mode (Fig. 1131 centre) with v — 0.33 mHz 
and rj = —0.95. The stability is due to the near extinction 
of the previous driving region and most of the significant 
energy contributing to damp the modes. The development 
and extinction of the region contributing to driving (with its 
maximum taking place at the location of the Z-bump) at the 
expense of a damping region is produced in the frequency 
range between ~ 0.2 to 1.5 mHz and is responsible for the 
oscillatory profile of the growth rate. 
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Figure 13. Left, Centre & Right: Differential worlc (solid line), running work integral (dotted line) and opacity (dashed line) for the 
gl90, gl 77 and p89 mode, respectively for a quadrupole mode of model 10. All plots have been scaled to arbitrary units. The vertical 
dashed-dotted lines in the right plot depict the convection zone. 
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Figure 11. Evolutionary track described by model 10 (initial 
diamond). Diamonds mark models found unstable and crosses 
those found stable. Error bars give the T ^ ff and logg dotermi- 
nation for J 1600+0748 by iFontaine et al.l feoOSh (bottom) and 
iRodrfeuez-Lopez et al.l ll2009a^ ftop). Extra mixing in the model 
carries helium to hotter regions were it burns giving rise to ther- 
mal pulses. The energy produced causes the star to expand and 
cool evolving to lower temperatures and when the burning finishes 
the star returns to higher temperatures. 

Finally, mode p89 (Fig. [I3l right) with ly = 10.03 mHz 
and Tj = —1.0 is representative of the high stability of low 
and high-radial order p-modes. The stability is due to a wide 
damping region at the location of the Z-bump. 

Model 10.1 of this same evolutionary sequence was 
found unstable in the frequency range (for I — 2, and simi- 
larly for other degrees); 

• 0.71 to 0.75 mHz (periods 22-23 min) which corresponds 
to modes from g81 to g76. 



6 MODEL 15: ETA 675T90MIXI 

This is the second model found unstable in our study. 
We found it drives two overstable modes. This driving is 
achieved through the action of three different driving re- 
gions: (i) the opacity bump at ~ 2, 000, 000 K which for 
this model is due to a combination of opacity from iron- 
peak elements and, because of their increased abundances 
due the extra-mixing, opacity from C and O, (ii) a driving 
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Figure 12. Growth rate vs. frequency for model 10. Diamonds in- 
dicate the modes for which the differential work and running work 
integral were plotted. The dashed-dotted lines delimit the two 
region of unstable modes. Higher frequency modes whore found 
highly stable and were not plotted for clarity reasons. 

region at ~ 2 10^ K, located immediately above the hy- 
drogen burning shell, and most likely due to the K-edge of 
highly ionised argon ([Rogers &i Iglesia3ll992l ). and (iii) an 
opacity bump that occurs at the transition from radiative to 
conductive energy transport. At the high temperatures and 
densities in the core the radiative opacity is dominated by 
free - free absorption. Because the core temperature gradient 
is relatively small, the radiative opacity increases with the 
density. On the other hand the electrons become more de- 
generate with increasing density and the conductive opacity 
decreases. The combination of these two opacity variations 
gives a maximum in the opacity, corresponding to the opac- 
ity bump in driving region (iii) (see Fig. llSp . Tests cancelling 
out the contribution of the nuclear energy generation show 
that the e mechanism is not responsible for this deepest driv- 
ing region, nor for region (ii). The instability was obtained 
for I = 2 modes (similar for other modes with other degrees) : 

• g27 with u = 2.08 mHz and = 0.07 and 

• g24 with u = 2.32 mHz and rj = 0.05. 

This corresponds to periods in the range 7-8 min. 
Model 15 has one of the lowest H mass fraction, and one 
of the largest He, C and O mass fractions (see Table [1]), 
which extend with constant values from the surface of the 
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Figure 15. Rosseland mean opacity ft (dotted line, left axis), and 
its opacity derivatives (right axis): kx (solid line) and ftp (dashed 
line). Red regions fulfil the criterion for instability. See text for 
details. 



star down in the envelope until about log q ~ —4. It has 
also the largest metallicity, Z=0.19, of the whole series of 
models. 

The propagation diagram and opacity of model 15 
(Fig. I14p shows that the region of mixed modes is reduced 
compared to that of model 8, which is representative of most 
of the models used in this study. We also note that the Z- 
bump has shifted to shallower regions, with its maximum 
at logg ~ —11 compared to \ogq ~ —10 for the mod- 
els with TcB = 70 000 K or logg ~ —8 for models with 
Teff = 45 000 K. 

A plot of the Rosseland mean opacity and its opacity 
derivatives, kt = 91nK/91nT|p and ftp = d\n K,/d\n p\t is 
shown in Fig. [15] Opacity derivatives play an important role 
in the k mechanism as t he necessary condi tion for driving is 
given by the expression l|Unno et al.lll989l ): 



Figure 18. Radiative (green) and conductive (blue) opacity vs. 
pressure. Temperature (maroon) is given in the right axis. The 
vertical dashed line marks the opacity maximum. See text for 
details. 



Regions in the star fulfilling this criterion are shown in red 
in Fig. 1151 All the regions involved in the effective driving of 
the two excited modes are potential driving regions, together 
with the region of the Z-bump, also involved in the p-mode 
excitation (see below) and the very surface of the star, which, 
on the contrary, has no significant weight in the driving. 

The age of this model from when it first became an sdO 
is about 10 000 yr, which is less than the linear growth rate 
e-folding times found: 7 10* and 10^ yr for ^24 and g27, re- 
spectively. However, the effective temperature remains near 
90 000 K for more than 10^ yr after this particular model. 
Therefore, from the rate of change of the effective temper- 
ature, the evolutionary time scale for this model is about 
10^ yrs, and modes may be able to develope observable am- 
plitudes. 

The normalised growth rate of the model is given in 
Fig. 1191 It has a low frequency region with 0.5 ^ v ^ 
3.0 mHz where the two overstable modes are driven, and 
in which high and low values of the growth rate alternate 
for consecutive modes. There is a highly stable region at 
intermediate frequencies 5.0 < < 13.0 mHz and finally, 
a tendency to driving is found in the high frequency range 
with 14.0 ^ 1^ ^ 20 mHz where low-radial order p-modes are 
found. 

An example of a quasi-unstable mode in the very low 
frequency region, below 1 mHz, is given in Fig. [16] (left) by 
the g63 mode with v = 0.93 mHz and 77 = —0.15. The high 
values of the growth rate are due to driving from region (i). 
At the high effective temperatures of this model this opacity 
bump occurs high in the envelope, at logg ~ —6, which 
favours pulsations, as heat interchange is more effective. 

An example of a mode in the low frequency region which 
has lost its tendency to driving is given in Fig. [16] (centre) by 
the g35 mode with 1/ = 1.63 mHz and ri = —0.72. In addition 
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Figure 16. Left, Centre & Right: Differential work (solid line), running work integral (dotted line) and opacity (dashed line) for the 
g63, g35 and g24 mode respeetively for a quadrupole mode of model 15. All plots have been scaled to arbitrary units. 




Figure 17. Left, & Right: Differential work integral (solid line), running work integral (dotted line) and opacity (dashed line) for the 
g6 and p7 mode respectively for a quadrupole mode of model 15. The vertical dashed-dotted lines in the right plot depict the convection 
zone. All plots have been scaled to arbitrary units. 



to driving region (i), a driving region at —1 < logq < —0.5 
(corresponding to ~ 8 10^ K in temperature) has developed 
near the location of the deep opacity bump due to the tran- 
sition from radiative to conductive energy transport in the 
semi-degenerate material. However, radiative damping right 
above this region exceeds the driving effect of both regions, 
damping the modes. 

One of the two unstable modes found, g24 with u = 
2.32 mHz and rj — 0.05 is shown in Fig. [TS] (right). A new 
driving region, which we attribute to Ar ionization, has de- 
veloped at logq ~ —3, where T ~ 210^" K. The model is able 
to achieve instability due to the combined driving of two of 
the three regions: the first and second, or the first and third 
regions, but modes are stable if only the first, or the joint 
second and third regions contribute to the differential work. 
It is the first time, to our knowledge, that driving of modes 
due to such a combined action of driving regions has been 
reported. This may give us a useful tool to explore physi- 
cal processes, in particular energy transport, in dense hot 
regions. 

A mode representative of the intermediate frequency 
region with the lowest values of the growth rate is the g6 
mode with ty = 6.51 mHz and rj = -1.0 (Fig. [TT] left). 
Modes at intermediate frequencies have maximum energy 
interchange at logg ~ —8.5, where the derivative of the 
opacity does not favour driving, which results in maximum 
stabilization of the modes. 

Fig. 1171 (right) shows mode p7with u — 18.05 mHz and 
r] = —0.45, as an example of a quasi-unstable mode in the 
high frequency region where the low-radial order p-modes 
are found. This mode falls within the region of observed 
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Figure 19. Growth rate vs. frequency for model 15. Diamonds 
indicate the modes for which the differential work and running 
work integral were plotted. 



modes for J1600+0748. Surprisingly, due to the high T^b of 
the model, the high values of the growth rate are due to a 
wide driving region caused by Z-bump. However, driving is 
not achieved as it is counterbalanced by a damping region 
located at the base of the Z-bump. 

Further evolved models in the evolutionary sequence 
were analysed and found stable. 
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r SUMMARY OF THE RESULTS 

We present a summary of the stability analysis carried out 
for all the analysed models and submodels in Table [2] We 
give Teff, logg, metallicity and if a tendency to instability 
- or actual instability - was found in the g- and/or p-mode 
spectrum and which driving region was responsible for it. 
In addition, we give the frequency and period range most 
favoured for instability. 

If we draw our attention first to the first batch of mod- 
els built for this study (models 1 to 9), all of them, except 
models 7 and 8, have a certain tendency to instability of the 
5^modes due to the gathering of driving energy at the loca- 
tion of the deep opacity bump at T ~ 2 10^ K. This opacity 
bump is mainly due to iron-peak elements in models which 
have not experienced mixing of the products of core helium 
burning into the envelope. Opacity from C and O plays a 
role for models in which mixing has occurred. Whether or 
not opacity from C and O is important for driving, we label 
this driving region C/O in Table (2] Modes more favoured 
for driving are high-radial order g^modes, i.e. long period 
modes, with typical frequencies of oscillation between ~0.5 
and 2 mHz, i.e. periods ~500 to 2000 s. We emphasis the 
unique behaviour of model 7 in which the tendency to in- 
stability of the (/-modes is caused by the Z-bump. 

We would expect that an increase in the C/O mass frac- 
tion, at least in the driving region, would lead to the eventual 
driving of the modes, as it is known that high metallicities 
help to increase the potential instability of the models. 

Models 1 and 8 have the Z-bump playing a role in the 
destabilization of the p-modefl Model 8 presented an oscil- 
latory profile of the growth rate with frequency, which lead 
us to analyse it in more detail, reveal ing mode trapping of 
g- and p-modes (|Rodn'guez-L6pez et al..,2009h ). 

While doing this study we were aware that standard 
sdB evolutionary models with solar and uniform metallic- 
ities were not able to drive modes, and c rude increases in 
the metallicity first (|Charpinet et al ]| 19961) and the so-called 
second generation models next ( Charpinet et al.] Il997bl ) 
were needed to excite pulsations and reproduce theoreti- 
cally the observed instability strip. The latter take diffusion 
(radiative levitation and gravitational settling) into account 
to produce a non-uniform profile of iron with depth, which 
is confirmed to be a key ingredient in driving pulsations 
jCharpinet et al.ll2009l '). 

Due to our negative results in exciting modes in this 
first batch of models, we aimed at building models as rich in 
metals as possible to test if it was possible to excite modes 
at some point via a k- mechanism associated with a C/O 
partial ionization zone. Therefore, a last batch of models 
(models 10 to 16) was built in which an ad hoc extra-mixing 
was added to increase the metallicity in the envelope. 

All these mixed-models have a tendency to driving of 
the high-radial order g^modes (see Table [2| , caused by the 
Z-bump (this was only found previously for model 7), the 
C/O partial ionization zone, or a combination of both, the 
C/O being responsible for the tendency to driving of modes 
with lowest frequencies. 



^ Models marked with a lowercase z in Table [2] mean tliat the 
model barely reached the 'tendency to instability', considered as 
such when the growth rate values of the modes reach at least -0.5. 



We finally achieved the instability of several models in 
the evolutionary sequence of model 10, that define a nar- 
row instability band in the range 45 000 ^ T^a ^ 54 000 K, 
4.2 ^ logg ^ 4.5, metallicity 0.14 < Z ^ 0.15 and periods 
around 20 min, 45 min and 1 h. The instability is caused 
by a classical K-mechanism due to the opacity bump of the 
iron-peak elements, that excite high-radial order ij-modes in 
the star. Further time-evolved models in the same sequence 
were found stable. 

Model 15, with T^s: = 90 000 K, log g = 5.86 and metal- 
licity Z=0.19, was found to have unstable high-radial or- 
der g^modes due to driving from a combination of a clas- 
sical K-mechanism in the C/O partial ionization zone at 
r ~ 2 10'^ K, a K mechanism operating at T ~ 2 10^ K 
probably caused by argon ionization, and a K-mechanism 
associated with an opacity bump that results from the tran- 
sition from radiative to conductive opacities. This model has 
one of the largest total metallicity, helium, carbon and oxy- 
gen mass fractions within the whole set of models. We found 
it drives two overstable modes in the period range 7-8 min. 

Models 10.5, 10.6 and 15, with temperatures between 
84 500 and 91000 K, show a tendency to instability of the 
low-radial order p-modes caused by the Z-bump. This be- 
haviour is similar to that found in model 1, also a high 
temperature model (Tgfr = 79 000 K) . This is a somewhat 
unexpected result, as for such high Tad the location of the 
Z-bump is quite high in the envelope of the star, where the 
driving efficiency is lower due to the lower density. For- 
tunately, this result is confirmed by actual driving of p- 
modes throug h a k mechanism ass ociated to the Z-bump 
is achieved bv lFontaine et al] (|2008l ) for two hot sdO mod- 
els with Tcff = 71 000 and 81 000 K, which include radia- 
tive levitation of iron. Those models have high tempera- 
tures fitting the first determination of physical parameters of 
J16004-0748, T cff = 71 000 ± 2725 K a nd logs = 5.93 ±0.11, 
also derived by iFontaine et al.l l|2008l ') . 



8 DISCUSSION AND CONCLUSIONS 

We carried out a non-adiabatic analysis of sdO equilibrium 
models to explore their feasibility as pulsators. We analysed 
27 sdO models out of 16 different evolutionary sequences 
and we found the first two sdO evolutionary models ever to 
drive high-radial order g-modes. In one model, with TaH = 
45 000 K, instability of g-modes is due to a classical k- 
mechanism in the partial ionization zone of the iron-peak 
elements at T ~ 2 10^ K. This shows that radiative levita- 
tion is not an absolute pre-requisite for pulsations of a model 
sdO star as long as the metallicity is increased to what are 
probably unrealistic values. Inclusion of radiative levitation 
of iron remains necessary for complete asteroseismological 
analysis of hot subdwarfs. In the second model, with TcS — 
90 000 K, g-modes are found to be unstable due to the combi- 
nation of driving from three distinct regions in the star. The 
outermost driving region is located at T ~ 2 lO*' K, where 
the iron-peak opacity bump is strongly augmented by con- 
tributions from ionization of carbon and oxygen, which we 
have have enhanced in abundance by ad-hoc extra-mixing 
during the evolution after the ZAHB. The second driving re- 
gion is located at T ~ 2 lO'^ K where a subtle opacity bump 
develops due to the ionization of argon. The innermost driv- 
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Table 2. Effective temperature, logarithm of surface gravity, current metallicity and region responsible for the tendency to driving — or 
actual driving — of the g- and/or p-modes spectrum. The last two columns give the frequency or period range of most favoured for 
driving (taking into account modes from Z = 1 to i = 4). The lower-case z for models 7 and 8 indicates the modes are under the limit of 
the tendency to driving. Models between horizontal lines have been calculated with small changes in the evolution code. Models in bold 
were found unstable. 



Model number 


Tefl 


log g 


Z 


g-modes 


p-modes 


Freq. range 


Period range 














(mHz) 


V') 


1 p675.8057 


79 000 


5.70 


0.07 


C/O 


Z 


0.5-1.5, 10-14 


700-2000, 70-100 


1.1 p675.9065 


90 000 


6.50 


0.07 


C/O 


- 


0.5-5 


200-2000 


2 p650.7960 


79 000 


5.95 


0.02 


C/O 


- 


0.5-1 


1000-2000 


3 etap 685 


55 000 


5.89 


0.02 


C/O 




0.5-2.5 


400-2000 


4 etap 690 


55 000 


5.95 


0.02 


C/O 




0.5-2.5 


400-2000 


5 etap 695 


55 200 


5.98 


0.02 


C/O 




0.5-3 


300-2000 


6 etap 700 


55 000 


6.02 


0.04 


C/O 




0.5-2.5 




7 eta 600t45 


45 000 


4.95 


0.05 


z 




2-5 


200-500 


7.1 eta600t70 


70 000 


5.95 


0.05 


C/O 




0.5-1.5 


700-2000 


8 eta650t45 


45 000 


5.26 


0.05 




z 


7-9.5 


100-150 


8.1 eta650t70 


70 000 


6.22 


0.05 


C/O 




0.5-3.5 


300-2000 


9 eta 700145 


45 000 


6.13 


0.07 


C/O 


- 


1-7 


150-1000 


9.1 eta 700170 


70 000 


6.34 


0.07 


C/O 




0.5-4 


250-2000 


10 eta 675t45mixi 


45 000 


4.18 


0.14 


z 




0.2-0.9 


1100-5000 


10.1 15773t54g45 


54 500 


4.54 


0.15 


z 




0.3-2 


500-3300 


10.2 15873t58g47 


58 300 


4.68 


0.14 


z 




0.3-2.5 


400-3000 


10.3 16973t69g51 


69 000 


5.14 


0.16 


c/o+z 




0.2-3.5 


300-5000 


10.4 17773t77g54 


77000 


5.39 


0.17 


C/0+ z 




0.2-1.5, 4-6.5 


700-5000, 150-250 


10.5 18273t84g56 


84 500 


5.63 


0.18 


C/O 


Z 


0.2-1, 8-11 


1000-5000, 90-125 


10.6 18573t91g59 


90 700 


5.93 


0.19 


C/O 


Z 


0.2-3, 17-21 


300-5000, 50-60 


11 ota 675t45mixnmi 


45 000 


4.36 


0.07 


z 




0.4-0.7 


1400-2500 


12 eta675t70mixli 


70 000 


5.16 


0.16 


C/O + z 




0.2-1, 2.5-4 


1000-5000, 250-400 


13 eta 675t70mix2i 


70 000 


5.20 


0.18 


C/O + z 




0.2-0.5, 3-4 


2000-5000, 250-300 


14 eta 675t70mixnmi 


70 000 


5.43 


0.07 


C/O + z 




0.4-0.8, 5-7 


1250-2500 


15 eta 675t90mixi 


90 000 


5.86 


0.19 


C/O 


Z 


0.3-3.5, 15-20 


300-3300, 50-65 


15.1 18742t98g63 


97 600 


6.28 


0.20 


C/O 




1.5-2 


500-700 


16 eta 675t90mixnmi 


90 000 


6.58 


0.07 


C/O 




2-4, 5.5-7 


250-500, 150-180 



ing region is located at the transition between radiative and 
conductive energy transport, and may be of interest as a tool 
to explore energy transport processes in dense hot regions. 

All of our sdO models, except one, present a certain 
tendency to instability for high-radial order (?-modes, with 
typical oscillation frequencies between ~0.5 and 2 mHz (i.e. 
periods ~500-2000 s). For most of the models this tendency 
to driving in usually due to the occurrence of a driving zone 
at the location of the C/O partial ionization zone. However, 
quasi-excitation due to the heavy elements partial ionization 
zone or a combination of both also occurred. 

The fact that we found actual driving due to a k- 
mechanism associated with C/O opacity in only one model 
may be due to the deep location of the C/O partial ioniza- 
tion zone within the star, where the high ratio of the thermal 
to dynamical time scales impedes efficient heat interchange. 
Hence, only hotter models with shallower partial ionization 
zones, such as model 15 for which we achieved actual driv- 
ing, may be able to overcome this problem. 

The tendency to instability of p-modes was not so com- 
mon, and was achieved for only 5 out of 27 models. We 
did not achieve actual driving of p-modes, those observed in 
the only sdO pulsator known to date, in any of our mod- 
els. In some of the hotter models, we found a tendency to 
driving of low-radial order p-modes, always triggered by the 



K-mechanism associated with the Z-bump. We partially at- 
tribute the inability to get actual p-mode driving to the shal- 
low location of the heavy elements partial ionization zone, 
where the density is low, making its contribution to the work 
integral not sufficient to drive the modes. Ho wever, this diffi- 
culty h as been now overcome by the results of lFontaine et al] 
(|2008l ) who have achieved the actual driving of p-modes for 
hot models which include radiative levitation and gravita- 
tional settling of iron. 

Finally, we wou ld li ke to remark that, as shown by 
iFontaine etal] (|2008l ') and lCharpinet et al.1 (|2009l '). inclusion 
of radiative levitation of iron in the models has revealed es- 
sential in achieving the actual excitat ion of p-modes. The 
determination of physical par ameters (|Fontaine et al.ir2008l : 
iRodn'guez-Lopez et al.ll2009a^ for the unique pulsating sdO 
star to date provides us with an excellent clue to pursue ob- 
servational searches for new sdO pulsators, with the aim of 
characterising their instability strip. Envisaged problems are 
the lack of catalogued sdOs with high effective temperatures, 
as historically, atmosphere grids reached only up to about 
55 000 K. The role of recent sky surveys such as the Sloan 
Digitally Sky SurvejQ and new tools as those offered by the 
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Virtual Observatorjjf] are crucial for spectroscopic and pho- 
tometric searches of new sdO candidates. The road ahead is 
challenging, as candidates will be mostly faint (magnitudes 
at least over 16 or 17) which will pose difficulties in achiev- 
ing high S/N light curves that allow to detect the lowest 
amplitude modes. Strategies as propose sdO candidates to 
be observed by asteroseismologic satellites such as CoRoT 
or KEPLER have to be considered. 
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